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Abstract 

The redshift-distance modulus relation, the Hubble Diagram, derived 
from Cosmological General Relativity has been extended to arbitrarily 
large redshifts. Numerical methods were employed and a density function 
was found that results in a valid solution of the field equations at all 
redshifts. The extension has been compared to 302 type la supernova 
data as well as to 69 Gamma-ray burst data. The latter however do not 
not truly represent a 'standard candle' as the derived distance modulii are 
not independent of the cosmology used. Nevertheless the analysis shows 
a good fit can be achieved without the need to assume the existence of 
dark matter. 

The Carmelian theory is also shown to describe a universe that is 
always spatially flat. This results from the underlying assumption of the 
energy density of a cosmological constant Q,a = 1, the result of vacuum 
energy. The curvature of the universe is described by a spacevelocity 
metric where the energy content of the curvature at any epoch is Qk ~ 
Qa — n = 1 — where il is the matter density of the universe. Hence the 
total density is always fix + — 1. 

Keywords: Cosmological General Relativity, high redshift type la super- 
novae. Gamma-ray burst, dark matter, distance modulus 



1 Introduction 

Carmeli's cosmology, also referred to as Cosmological General Relativity (CGR), 
is a space-velocity theory of the expanding universe. It is a description of the 
universe at a particular fixed epoch of cosmic time t and involves a new dimen- 
sion V, the velocity of the expansion of the fabric of space itself, in which the 
galaxies are sitting. For a full introduction to the theory, the reader should 
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become familiar with either of the late Moshe Carmeli's books Cosmological 
Special Relativity, 2nd Edition [2] or Cosmological Relativity [4] . 

In the 1920s Edwin Hubble discovered that the observed redshift z in the 
light emitted from a distant source is directly proportional to the distance D to 
the source, viz. v = Hq D, where Hq is Hubble's constant. Carmeli incorporated 
this law as a fundamental axiom, with Hq^ « r, a universal time constant, 
into a general 4D Riemannian geometrical theory satisfying the Einstein field 
equations (Ref. appendix A). The observables, redshift and distance (Hubble 
distance) , are the coordinates of his theory. 

The expanding universe is then described by its spacevelocity "phase space" 
coordinates {v,r,9,(j)), which do not involve time t but this new dimension v. 
Cosmic time is measured from the present epoch, with t ~ Q, back toward the 
beginning with t — t, but determined from the expansion velocity of the fabric 
of space, i.e t — > r as w — > c. In this case the galaxies act as tracers embedded 
in the expanding space. This is different to the approach in the 4D Friedmann 
theory with coordinates {t, r, 6, <j)). But Carmeli never solved the full 5D theory 
involving the time dimension as well. In section[S]I show how the time dimension 
is included in a general 5D CGR theory for the universe. Surprisingly, or may 
be not, what results is that the universe is Euclidean at all epochs. 

Usually the Hubble law is stated in terms of the velocities of the receding 
galaxies, when actually those velocities are determined from the measured red- 
shifts of spectral lines in the galaxy light. Carmeli used this fact to formulate his 
new theory. In practice redshift, and not velocity, is what astronomers actually 
measure, as well as distance determined from the luminosity of the source. In 
section [3] I describe how the CGR theory can be compared with observations 
using luminosity distance as an independent variable, and how this is different 
to FRW theory. In section H] I discuss how to extend the theory to a much 
greater range of validity in redshift space by solving Carmeli's equation directly 
(though numerically) for redshifts, instead of working in velocity space. In sec- 
tion O I show how well the theory fits the observations, including both type la 
supernova and Gamma-ray bursts. 

2 Phase space equation 

Carmeli initially developed his special theory from considerations of a universe 
without matter, then later solved the field equations for a fully covariant general 
4D theory. He noticed the analogue with special relativity. In Cosmological 
Special Relativity (see chap. 2 of [2]), instead of the Minkowski line element of 
special relativity we can write the line element as 

ds^ ^T^dv^ - dr^, (1) 

where dr^ = dx^ -\- dy^ + dz^. Now by equating ds = it follows from ^ that 
rdv = dr assuming the positive sign for an expanding universe. This is then the 
Hubble law in the small v limit, where r is the Hubble-Carmeli time constant, 
the inverse of the Hubble constant in the limit of weak gravity. 
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One sees the similarity, where instead of the hght cone of special relativity, 
we get the concept of a velocity cone. Galaxies are distributed in the universe 
such that as their recession speed approaches the speed of the expansion of space 
they approach the limits of the cone. So, in general, this theory requires that 
the expansion of the universe be described by ds = 0. The following describes 
the result when matter is considered in the general theory. 

The line element in CGR [5] 

ds^ = T^dv^ - e^dr^ - R^{d0^ + sin^Odcj)'^), (2) 

represents a spherically symmetrical isotropic universe, that is not necessarily 
homogeneous. 

Using spherical coordinates (r, 0, </>) and the isotropy condition d6 = dcj) = 
in ([2]) then dr represents the radial co-ordinate distance to the source and it 
follows from ([2]) that 

T^dv^ ~ e^dr^ = 0, (3) 
where ^ is a function of v and r alone. This results in 

where the positive sign has been chosen for an expanding universe. 

Carmeli found a solution to his field equations, modified from Einstein's, 
(see [8] and [21 [5] ) which is of the form 

with R = r and hence R' — 1, which must be positive. From the field equations 
and ([5]) we get a differential equation 

/' + ^ = -i^T'^Peffr, (6) 

where /(r) is function of r and satisfies the condition f{r) + 1 > 0. The prime is 
the derivative with respect to r. Here k — SttG/c^t'^ and pe// = p — Pc where p 
is the averaged matter density of the Universe and pc = S/SttGt^ is the critical 
density. Finally the solution to the inhomogeneous equation ^ was found to 
be 

/(r) = -^rW' = (l-f^)^- (7) 
Therefore this requires that 

In the above Carmelian theory it initially assumed that the Universe has 
expanded over time and at any given epoch it has an averaged density p, and 
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hence Peff- The solution of the field equations has been sought on this basis. 
However because the Carmeli metric is solved in an instant of time (on a cos- 
mological scale) any time dependence is neglected. In fact, the general time 
dependent solution has not yet been found. But since we observe the expanding 
Universe with the coordinates of Hubble at each epoch (or redshift z) we see 
the Universe with a different density p{z) and an effective density peff{z). In 
section |4] this fact is taken into account to extend the validity of the solution. 




3 Comparison with observation 

In order to compare CGR theory with observations one requires a relationship 
that gives proper distance as a function of redshift. From ([4]), ([5]) and ([7]) we 
get the differential equation 



(9) 

where il = p/ Pc- Equation ([9]) relates the velocity of the receding sources to their 
distance. Carmeli integrated ([9]) and found the redshift distance relationship 

r _ smh ipVT^) 

CT ' ^ ' 

where, 

V {l + zf-l 
c (1 + z)"^ + 1 

and the average matter density of the universe, a function of redshift z, must 
be found. 

This equation was used in both Hartnett [S], which forms the basis of the 
work presented here, and Oliveira and Hartnett [13] where comparisons were 
made to the high redshift SNe la data from Riess et al [16] , Astier et al |lj and 
Knop et al [12]. The proper distance is converted to magnitude as follows. 

m{z)^M + 5log[DL{z;n)], (12) 

where is the dimensionless "Hubble constant free" luminosity distance. Refer 
[TUdl]. Here 

M^5log(^-^— ^ +25 + Mb + a. (13) 

The units of cr are Mpc. The constant 25 results from the luminosity distance 
expressed in Mpc. However, M in represents a scale offset for the distance 
modulus (m-Ms). It is sufficient to treat it as a single constant chosen from the 
fit. In practice we use a, a small free parameter, to optimize the fits. 

In CGR the luminosity distance was found [10] to be slightly different to the 
expression used in the FRW theory. This fact, results from a relativistic effect 
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on the emission times of the photons from the distant source. In CGR, times 
at cosmological distances add according to a relativistic addition law, [3] when 
referred to the observer a,t t — 0. When this is taken into account we get an 
additional factor (1 — /3^)~^/^, hence Vl is given by 

Pi + (1-/32)-^/' (14) 

using (|10p . which is a function of f2 and z. 

The parameter M incorporates the various parameters that are independent 
of the redshift, z. The parameter Mg is the absolute magnitude of the supernova 
at the peak of its light-curve and the parameter a allows for any uncompensated 
extinction or offset in the mean of absolute magnitudes or an arbitrary zero 
point. It also acts as a correction to the value of t since M contains it. The 
absolute magnitude then acts as a "standard candle" from which the luminosity 
and hence distance can be estimated. 

The value of Mb need not be known, neither any other component in M, 
as Ai has the effect of merely shifting the fit curve ([U]) along the magnitude 
axis. However by choosing the value of the Hubble-Carmeli constant r = 4.28 x 
lO'^^ s = 13.58 Gyr, which is the reciprocal of the chosen value of the Hubble 
constant in the gravity free limit h — 72.17 ± 0.84 (statistical) km.s~^Mpc~^ 
^M= 43.09 + Mb + a. 



4 Extended redshift range 



The proper distance from (jlOp however has limited application because it has 
essentially been found by integrating ^ assuming constant density then using 
(jlip . The result being only valid over a limited redshift range. In order to 
improve on this, because the observables are distance (magnitude) and redshift 
not velocity, the differential equation ([9]) must first be converted to dr/dz using 
the chain rule as follows. 



1 f dr dz\ 1 ~ fl{z) 2 



r^\dzdv) =1 + ^^^^' (15) 

where it must be remembered that matter density is also a function of redshift. 
The derivative dv/dz is then calculated from ([TI]) and substituted into 
which becomes 

1 fdrV _ 16(1 + z)2 



dzj (l + (l + z)2)4 



l + {l-niz)) . (16) 



This differential equation is then solved for r/cr, but because it involves mixed 
terms on the rhs, it must be numerically solved. But we can employ a 'boot- 
strapping' technique. By substituting the hmited redshift solution from (fTU)) 
into the rhs and assuming the averaged matter density in the universe 

n{z)^n,nii + zf, (17) 
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equation ([TC]) was numerically solved for r/cr. The result is shown in curve 1 
of fig. [TJ This is compared with r/cr resulting from the initial equation (|10p 
represented by curve 2. 

Then r/cr (of curve 1) was then taken as input on the rhs of and 
iterated again to numerically solve for r/cr. However because of the requirement 
^ on the solution from the field equations the result is only valid to about 
z = 2.5 with the assumed density dependence (jl7p on the present epoch density 
= 0.04. Dimensionally the matter density depends on the length cubed if 
we assume constant mass, but the precise density dependence is unknown at 
high redshifts. Therefore curve 1 of fig. [1] plotted beyond z = 2.5 can only be 
considered approximate. 

Now in CGR there is the scale radius 

Ro = ^^=- (18) 

From and the condition ([5]) the range of validity of the solution © is 
< r/cr < 1 < Rf) if Q < 1 and < r/cr < i?o if > 1. This means the 
solution is valid for all values of the matter density ^{z) < 2. 

Now if we again use the solution (curve 1 of fig. [T]) and solve p6|) for matter 
density fl, it results in a density function such that ([8]) is always true. The 
density function is shown as curve 1 in fig. [2] and is compared with the initial 
density function (fT7|l shown as curve 2. Finally using this as the density function 
in (|16p and numerically solving for r/cr results in curve 3 in fig. [1] 

The new density function has some unexpected features, but it produces a 
smooth monotonically increasing function of distance on redshift as expected. 
See fig. [3] where it has been plotted to z = 100 for the case where ft,n = 0.04. 
The fluctuating part of il{z) for z > 4 is approximately 



1 + 1.1 cos^ 



/3v/0.04(l + z)3- 1 + 1.6 , (19) 



where /3 is determined from (jlip . The function (|19p has a decreasing period 
as a function of a minimum a,t — 1 and a maximum slightly greater than 
fl = 2. Any density function that remains within the range of this function for 
z > 4, even a constant density, will yield a smooth monotonic redshift-distance 
relation. However the oscillating density as a function of redshift could represent 
episodic creation from the vacuum as the universe expands. Particle production 
at those epochs is a byproduct of the expansion process. Gemelli [TT] found 
from a hydrodynamic solution of the 5D problem that particle production must 
occur from the expansion provided the Universe is not isentropic. 



5 Quality of curve fits 

The initial equation (fTO|) has been curve fitted to SNe la data [l3j as described 
in section[3]and excellent fits resulted without any dark matter component. The 
best statistical fit resulted in a value of flm — 0.0401±0.0199, which is consistent 
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with the observed baryonic matter density [7]- The fit was accompHshed with 
the combined Gold and Silver SNe la data sets of Riess et al [16] with that of 
Astier et al [T] that extended to z = 1.75, the current limit of observational 
data. Inspection of the r/cr curve (curves 1 or 3 shown in fig. [T]) and the initial 
curve (curve 2) for the region 1.5 < z < 1.75 indicates that only the last two 
SNe la data would be sensitive to the numerically determined extension here. 
See crosses in fig. |4l 

Using r/cT from the extended equation I have curve fitted to the data 
set of Riess et al and also to combined data set of Riess et al and Astier et 
al. However before using Astier et al data I added 0.14 magnitudes to their 
quoted errors to align their errors with those of Riess et al. Otherwise Astier et 
al seem to have underestimated their magnitude errors, which unfairly bias the 

calculations. The resulting y^s for the model here are shown in Table I. 

Clearly more SN la data, at higher redshifts, are needed to test this predic- 
tion. Consequently GRBs have been used to get distance modulus and redshift 
as independent parameters (from Schaefer [T^) and are shown as dots in fig. 
21 However the magnitude calculation is not really independent of the chosen 
cosmology as explained in the latter reference. Nevertheless using the data from 
columns 2 and 8 of Table 6 TT , and making the necessary changes to convert the 
magnitudes from the FRW Concordance Model to that of the Carmeli model, 
new distance modulii are obtained. 

In order to use the data from Schaefer [T7] first the specific details particular 
to a cosmological model he used are subtracted then the specific details from 
the CGR theory are applied. This means small corrections are made and are 
shown in fig. [6] as a function of redshift. This also means the magnitudes 
calculated this way are not really independent of our chosen CGR cosmology 
either. In fact, problems related to the cosmological dependence of the GRBs 
calibration, in principle, may invalidate such a procedure. Nevertheless we may 
get some benefit from this exercise. In figs |4] and [5] both SN la and GRB 
data are represented on the same plot. However the SN la magnitudes are not 
dependent on any particular cosmology and present a much stronger constraint 
on the applied theory. Because it is shape of the curve fitting to the data that 
carries the most information if the theory fits the SN la data it should also fit 
the GRB data. 

The following are the corrections applied to convert the distance modulus 
of Schaefer to that of the extended theory of CGR developed here. Since the 
luminosity distance in CGR is different to that in FRW theory, by the factor 
(1 — it's effect must be added to the distance modulus. Therefore the 

luminosity distance equation (10) of Schaefer must be replaced by that in CGR. 
So equation (|14p . with r/cr determined from the numerical solution of (jl6p and 
— 0.04, was used to calculate a correction, which was added as magnitude 
to the distance modulus. Equation (10) of Schaefer [T7] with Q.^ — 0.27 and 
fiA = 0.73 for the Concordance Model was numerically integrated and removed 
as a magnitude from the data. These three corrections have the combined effect 
on the data as shown in fig. [6l 

Finally by fitting the extended equation (fT6|) with ilm = 0.04 to this data 
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and to the combined set of Riess et al and Astier et al it was determined that 
0.1988 magnitudes needed to be added to shift the GRB set so it could be 
represented on the same plot as the SN la data. This was done by determining 
the best fit value of a from ([TH)) for each data set then adding the difference to 
the GRB distance modulii so a is the same for both. Hence a fit to the total 
data set can then be accomplished. See fig. There the data and fits using the 
extended Carmeli theory with = 0.02, 0.04 and 0.27 are shown on a redshift 
log axis. This better shows the low redshift fit. At high redshifts the fits with 
the smaller values of f2m are favored. 

In order quantify the goodness of the least squares fitting I have used the 
parameter which measures the goodness of the fit between the data and 
the theoretical curve assuming the two fit parameters a and Q.m- Hence is 
calculated from 

^ 1 

X" = ^ — [(m - M){z), - (m - M)(z„fc,),]' , (20) 
1=1 

where N are the number of data; (m — M){z) are determined from p2p with 
fit values of a and (™ — M)[zobs) are the observed distance modulus data 
at measured redshifts Zofcs! are the published magnitude errors. The values 
of /N (w Xd o f) ^^"^ shown in Table I, calculated using published errors on 
the distance modulus data. The best fit value of a is found at low redshift and 
then applied to all fits with various values of 



Table I: Curve fit parameter 



Data set(s) 




XVN 


SNe la 


0.00 


1.4120 


Riess et al Gold & Silver 


0.02 


1.4047 


N = 185 


0.04 


1.3987 


a = 0.2400 


0.06 


1.3941 


(a from z < 0.23) 


0.10 


1.3888 




0.27 


1.4258 


SNe la 


0.00 


1.2366 


Riess et al + Astier et al SNLS 


0.02 


1.2341 


N = 302 


0.04 


1.2326 


a = 0.1894 


0.06 


1.2321 


(a from z < 0.24) 


0.10 


1.2340 




0.27 


1.2871 


SNe la + GRBs 


0.00 


1.2191 


Riess + Astier + Schaefer 


0.02 


1.2005 


iV = 371 


0.04 


1.1902 


a = 0.2113 


0.06 


1.1860 


(a from z < 2) 


0.10 


1.1870 




0.27 


1.2538 



From Table I the smallest x^/N value are generally obtained where 0.04 < 
r2„j < 0.06, especially when the GRBs are included. But the quality of the 
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fits are not strongly dependent on the matter density. The locally measured 
baryonic matter budget has determined that 0.007 < SI;, < 0.041 [7] where a 
Hubble constant of 70 km.s^^ Mpc^^ was assumed. This region is consistent 
with the most probable values of fim from this analysis. Definitely a matter 
density of ^Im = 0.27 is not necessary for the best fit. Therefore no exotic dark 
matter need be assumed. However to properly test this model much more SNe 
la data are needed where z > 2 and much better constrained magnitude errors. 
The GRB data are not determined independently of the cosmology tested and 
therefore by themselves don't provide a 'standard candle' test of the theory. 



6 Spatially flat universe 

In section[2]thc time part of the metric g^i^ was ignored. Here we use coordinates 
a;^ — (ct, r, 6, cf), tv) where ^ = — 4, but the full 5D description of cosmos 
has never been determined in the Carmelian theory. So the question must be 
answered, "What is the goo metric component for the large scale structure of 
the Universe in CGR?" 

First note from ([5]) and ([6]) the gn metric component 

511 = - fl + V . (21) 



Using the scale radius we can define an energy density from the curvature 



which, when we use ([T5|. becomes 

nK = l-fl. (23) 

This quantifies the energy in the curved spacevelocity. 

In the FRW theory the energy density of the cosmological constant is defined 
PA = A/SttG hence 

A 



- (24) 

Even though the cosmological constant is not explicitly used in CGR, it follows 
from the definition of the critical density that 

' ' (25) 



when the cosmological constant A is identified with Therefore in CGR it 

follows that 

A 

3^2 



f^A^TTTI -A(-) =1. (26) 
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In CGR h = T~^. This means that in CGR the vacuum energy pyac = A/SttG 
is encoded in the metric via the critical density since pe/ f = P ~ Pc principally 
defines the physics. So JIa = 1 identically and at all epochs of time. (The 
determination of IIa in [S] was flawed due to an incorrect definition.) Also we 
can relate f^A to the curvature density by 

nK = nA-n, (27) 

which becomes 

nk = nA-nrr,, (28) 

at the present epoch (2 w 0). Here fl ~ i7„j(l + z)'' {z ^ 1) and hence fl^ ilk 
as z 0. 

Finally we can write for the total energy density, the sum of the matter 
density and the curvature density, 

fit ==n + nK ^^ + 1-^ = 1, (29) 
which means the present epoch value is trivially 

This means that the 3D spatial part of the Universe is always flat as it ex- 
pands. This explains why we live in a universe that we observe to be identically 
geometrically spatially fiat. The curvature is due to the velocity dimension. 
Only at some past epoch, in a radiation dominated universe, would the total 
mass/energy density depart from unity. 
I therefore write 

<?oo(r) = l + (l-r!t)r2, (31) 

where r is expressed in units of cr. Equation (|3ip follows from goo = 1 ^ 4$/c^ 
where $ is taken from the gravitational potential but with effective density, 
which in turn involves the total energy density because we are now considering 
spacetime. The factor 4 result from a comparison of the CGR theory with 
Newtonian theory. 

Clearly from (|29|) it follows that 500(f) = 1 regardless of epoch. Thus from 
the usual relativistic expression 



1 + Zgrav = \ I ~ 1, (32) 

V 9oo[r) 

and the gravitational redshift is zero regardless of epoch. As expected if the 
emission and reception of a photon both occur in fiat space then we'd expect 
no gravitational effects. 

Therefore we can write down the full 5D line element for CGR in any dynamic 
spherically symmetrical isotropic universe, 

ds^ = c^df -(1+ ^^-T^rA dr^ + T^dv^. (33) 
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The 9 and (p coordinates do not appear due to the isotropy condition d9 = d(p = 
0. Due to the Hubble law the 2nd and 3rd terms sum to zero leaving dt — ds/c, 
the proper time. Clocks, co-moving with the galaxies in the Hubble expansion, 
would measure the same proper time. 

However inside the Galaxy we expect the matter density to be much higher 
than critical, ie ilgaiaxy ^ 1 and the total mass/energy density can be written 

^o\galaxy — ^galaxy ~ ^galaxy : C^"^) 

because fife « 1, since it is cosmologically determined. Therefore this explains 
why the galaxy matter density only is appropriate when considering the Poisson 
equation for galaxies. [S] 

As a result inside a galaxy we can write 

900{r) = I + ^K-^J^ + ^galaxy-^J^, (35) 

in terms of densities at some past epoch. Depending on the mass density of the 
galaxy, or cluster of galaxies, the value of goo here changes. As we approach 
larger and larger structures it mass density approaches that of the Universe as 
a whole and 500 ^ 1 as we approach the largest scales of the Universe. Galaxies 
in the cosmos then act only as local perturbations but have no effect on il.K- 
That depends only on the average mass density of the whole Universe, which 
depends on epoch (z). 

Equation psp is in essence the same expression used on page 173 of Carmeli 

in his gravitational redshift formula rewritten here as 



A2^ l l + nKrl/c^T^-Rs/r2 , . 

Ai ^ l + QKrl/c^T^ -Rs/ri- ^ ' 

involving a cosmological contribution {Hkt^ /c^t^) and Rs = 2GM/(? , a local 
contribution where the mass M is that of a compact object. The curvature 
{^k) results from the averaged mass/energy density of the whole cosmos, which 
determines how the galaxies 'move' but motions of particles within galaxies is 
dominated by the mass of the galaxy and the masses of the compact objects 
within. Therefore when considering the gravitational redshifts due to compact 
objects we can neglect any cosmological effects, only the usual Schwarzschild 
radius of the object need be considered. The cosmological contributions in (I36p 
are generally negligible. This then leads back to the realm of general relativity. 



7 Conclusion 

The solution of the Einstein field equation using the spacevelocity of Moshe 
Carmeli requires the condition ([5]). This means that the redshift-distance rela- 
tion can only be valid, in general, where O < 2, in a matter dominated universe. 
If follows then if the Universe underwent periodic particle creation at high red- 
shift, while the expansion progressed, then we can construct a Hubble Diagram 



11 



that is valid over arbitrary redshift. This approach extends the Carmeli theory 
to where it makes a clear prediction and can be tested with future space-based 
telescopes, like the James Webb. 

The current plethora of cosmologies seem to require a few extra degrees of 
freedom, for example, scalar fields, dark energy or dark matter, but are they 
really capable of being distinguished one from another when tested against ob- 
servations? If a future space-based telescope was able to resolve well enough type 
la supernovae to redshifts approaching z — 10, then these various cosmologies, 
which fit the same data sets at much smaller redshifts, may be distinguished 
one from another. I believe only at these very high redshifts, provided the un- 
certainties on the data are small enough, can we hope to accomplish this task 
and really test Carmeli's cosmology. 

However using the Carmelian metric we are able to answer the question, 
'Why do we observe the Universe at this present epoch to be Euclidean?' 
Carmeli's choice of metric implicitly involves a vacuum energy density, such that 
when the matter density of the Universe reaches that of the vacuum density the 
fabric of space is relaxed and the expansion coasts (i.e. is unaccelerated) . The 
requirement on the energy content of the vacuum leads to the Universe being 
spatially flat at all epochs. This is because the total curvature of the Universe 
{flk = 1 — i^m at the present epoch) is the result of the average matter content 
warping spacevelocity away from flatness whereas the 3D spatial part remains 
flat. So at the present epoch Jlo = + = 1- 
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Figure 1: Redshift distance relation: Curve 1 represents r/cr from the solution 
of (fT6|) with ri(z) = r2m(l + z)^, curve 2 represents r/cr from (fTO|) with = 
rj„j(l + z)-^ and curve 3 represents r/cr from the solution of ^TE\i but with the 
density taken from curve 1 of fig. [2l 
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Figure 3: Redshift distance relation to z = 100 using density function of curve 
1 from fig. [2l 
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Figure 4: Distance modulus data from the combined SN la data set of Riess et 
al and Astier et al (crosses) and the GRB data (dots) derived from Schaefer [17] 
with curve fits using the extended equation ()16|) . The top solid curve represents 
the distance modulus where fim — 0.02, the broken curve below it is where 
rim = 0.04 and the dashed curve below the latter is where ilm = 0.27. 
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Figure 5: All 371 distance modulus data from the combined SN la data set of 
Riess et al and Astier et al and the GRB data derived from Schaefer Tf] with 
curve fits as shown in fig. |4l but on a redshift log scale. The top solid curve 
represents the distance modulus where 17„i — 0.02, the broken curve below it is 
where J7m = 0.04 and the dashed curve below the latter is where = 0.27. 
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Figure 6: The correction applied to the distance modulus data of Schaefer [IT] 
to account for the Carmeli cosmology. 
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